Spin currents in rough graphene nanoribbons: 
Universal fluctuations and spin injection 
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We investigate spin conductance in zigzag graphene nanoribbons and propose a spin injection 
mechanism based only on graphitic nanostructures. We find that nanoribbons with atomically 
straight, symmetric edges show zero spin conductance, but nonzero spin Hall conductance. Only 
nanoribbons with asymmetrically shaped edges give rise to a finite spin conductance and can be used 
for spin injection into graphene. Furthermore, nanoribbons with rough edges exhibit mesoscopic 
spin conductance fluctuations with a universal value of rmsGs ~ 0.4e/47r. 



PACS numbers: 85.75.-d 73.63.-b 72.25.-b 73.22.-f 

After their experimental discovery in 2004 [1 , mono- 
layers of graphite have attracted much experimental and 
theoretical attention owing to their unusual band struc- 
ture [2^ . Graphene has also been suggested as a good can- 
didate for spin based quantum computing and spintron- 
ics [3], as it is expected to have long spin decoherence/re- 
laxation times 0]. This prospect led to the recent 
interest in generating and manipulating net spin dis- 
tributions in graphene. Recently, spin injection from 
ferromagnetic metal contacts into graphene has been 
achieved El H E] • 

Transport properties of graphene nanoribbons (GNR) 
are expected to depend strongly on whether they have 
an armchair or zigzag edge [9 . In GNRs with zigzag 
edges, transport is dominated by edge states which have 
been observed in scanning tunneling microscopy p!0^. 
Moreover, owing to their high degeneracy, these states 
are expected to be spin polarized [11], making zigzag 
GNRs attractive for spintronics [12]. In addition, edge 
states are expected to occur also in nanoribbons with 
other edge orientations [13] . Recently, the first transport 
experiments have been performed in narrow ribbons of 
graphene [14], albeit with not well defined edges. Recent 
theoretical work focused on charge transport through 
rough GNRs [TF , but spin transport properties have not 
been explored yet. 

In the present work, we focus on spin transport in 
GNRs with rough zigzag edges. Ideal zigzag GNRs are 
not efficient spin injectors due to the symmetry between 
the edges with opposite magnetization. In order to obtain 
net spin injection, this symmetry must be broken. Exist- 
ing proposals to achieve this require very large transverse 
electric fields [12 . We sidestep this difficulty by showing 
that edge imperfections (such as vacancies), which usu- 
ally cannot be avoided experimentally, break the symme- 
try between the edges and lead to a finite spin conduc- 
tance of the GNR. Thus, rough zigzag GNRs can be used 
as spin injectors or detectors in graphene spintronics. 

We start with a description of the electronic ground 
state properties of the zigzag GNR, which captures the 




FIG. 1: (color online) Ground state spin density for (a) an 
ideal and (b) an imperfect zigzag GNR. Blue (red) corre- 
sponds to up (down) spin density, (c) Band structures of an 
ideal GNR obtained from DFT and tight-binding approaches. 

essential physics relevant to spin transport, given by the 
single band tight-binding Hamiltonian [TT] 

H = Yl^ij4,s^j,s -^^^i- 4,s^s,s'C^^s' • (1) 
ij,s i,s,s' 

Here tij = t if i and j are nearest neighbors, tij = if i 
and j are next nearest neighbors [16 , and a are the Pauli 
matrices corresponding to the spin degree of freedom. 
The local magnetization can be obtained from the 
self consistency condition or ab initio calculations. 

Our ab initio results, obtained using the spin-polarized 
density functional formalism (DFT) [17 , agree with the 
reported finding [TTJ |l2l [18] that the local magnetization 
is staggered in the electronic ground state, as shown in 
Fig. [TJa). At zero doping the antiferromagnetic (AF) or- 
dering generates a gap in the single particle spectrum. 
We now dope the GNR in order to move into a regime 
with open conduction channels. This can be achieved in 
practice by applying a gate voltage or chemical doping. 
Our DFT results indicate that a finite amount of dop- 
ing reduces the AF gap and the local magnetization, but 
does not destroy the AF ordering. We obtain the critical 
value of this doping as ^0.5 electrons {^OA holes) per 
zigzag edge atom. Furthermore, our DFT calculations 
show that not only perfect, but also rough zigzag rib- 
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FIG. 2: (color online) Spin injection profile from (a) an ideal 
GNR and (b) a GNR with a distorted edge into a region of n- 
doped graphene. Nonequilibrium densities for spin up (down) 
electrons are shown in blue (red). 

bons exhibit spin polarization (Fig. [ijb)). In addition, 
the formation of multiple spin domains at zigzag edges 
is energetically prohibitive. In summary, our DFT cal- 
culations show that it is possible (i) to dope the GNR to 
make them conductive and (ii) to introduce disorder at 
the edges while retaining the magnetic ordering. 

Next, we further simplify the mean field description of 
Eq. ([T]) by ignoring the variation of within a sublat- 
tice. A spatial dependence of changes the amount of 
band dispersion, modifying the energy window, within 
which the transport predominantly involves the edge 
states. This leads to the single particle Hamiltonian 

H,nf = e{k)Ti + A{k)T2 + A{k)l + m • ars, (2) 

where e(/c), A(/c) and A{k) are obtained by Fourier trans- 
forming Eq. (1), and are the Pauli matrices correspond- 
ing to pseudospin(sublattice) degrees of freedom [19]. 
The AF exchange field m is obtained by fitting the band 
structure to DFT results (see Fig. [TJc)). 

In the following, we focus on transport properties of the 
GNR. We work in the linear response regime so that all 
the transport properties of the GNR are specified by the 
effective single-particle Hamiltonian ([2|. The spin con- 
ductance [20^ of a GNR is given by Gs = (e/47r)(T| -T|), 
where is the transmission probability for spin up 

(down). The conducting channels with energies clos- 
est to the Fermi energy of the undoped system reside 
on a single sublattice and are fully spin polarized ow- 
ing to the staggered magnetization. These states are 
extended along the ribbon axis, but localized near the 
(zigzag) edges, with the spin up channel localized at one 
edge and the down channel on the opposite edge. The 
transverse localization length of these states depends on 
their Fermi momentum kp that may be modified by shift- 
ing the Fermi energy Ep. As one moves away from the 
X point, the transverse localization length increases as 
-^edge ^ — a/ ln(2 cos(A:Fa/2)), where a = 2.46 A is the 
hexagonal lattice constant ^ . Owing to the spatial sep- 
aration of the edge states, the scattering of spin up and 
spin down carriers occurs only at the edge, where they 
reside, and is unaffected by the opposite edge. Distin- 
guishing a left (1) and right (r) edge of the nanoribbon, we 



approximate T|(|) by Ti(r), where Ti(r) is the transmission 
probability of the corresponding edge state, assuming the 
opposite edge is not disordered. The transport properties 
of the zigzag GNR are thus essentially those of two inde- 
pendent wires, oppositely spin polarized and connected 
in parallel between the reservoirs. We note that previous 
studies of edge state transport ^ assumed vanishing 
next nearest neighbor hopping t', and obtained results in 
apparent contradiction to the picture presented above: 
If t' were zero, the charge density would be localized at 
the edges, but the current density would be extended 
through the GNR. This leads to the incorrect conclusion 
that edge states would scatter equally from impurities 
at both edges. In reality, the edge states show non-zero 
dispersion (such as due to t' 0). In this case, the cur- 
rent flow is also localized at the edges [22J validating the 
two-wire model, as we show below. 

For an ideal, impurity-free GNR, we have Ti = T^, 
which leads to vanishing spin conductance. This is con- 
firmed by quantum transport simulations [23 and an il- 
lustrative example is shown in Fig.|2|a): Both edge chan- 
nels transmit equally. However, as the edge states enter 
the bulk graphene, they are deflected: In the GNR, the 
pseudospin is predominantly in z-direction and tied to 
the electron spin, whereas in the bulk pseudospin is in- 
plane and tied to the current direction. At the interface, 
the 2;-component splits into states with positive and neg- 
ative velocity perpendicular to the boundary. The state 
with velocity towards the boundary is scattered [22] and 
thus, upon entry, states at opposite edges (which carry 
opposite spins) deflect in opposite directions, leading to a 
finite spin Hall conductance (Fig.[2|a)). Finite spin con- 
ductance can be obtained, however, for imperfect GNRs: 
An obstacle scatters the spin channel localized at the 
same edge more effectively, leading to a non-vanishing 
spin conductance and spin injection (Fig.|2|b)). Whereas 
the efficiency of the spin Hall effect is limited by the mean 
free path, and thus ballistic microstructures are needed 
to observe it, the efficiency of spin injection with edge 
defects is limited only by the spin relaxation length and 
can be used to inject spins into diffusive systems. 

From an experimental perspective, unless the GNRs 
are specifically fabricated with edges of different rough- 
ness, the average conductance of both spin channels is 
equal, quenching the ensemble-averaged spin conduc- 
tance. Yet, in the mesoscopic regime, sample-to-sample 
fluctuations of ^ lead to a non- vanishing variance of 
the spin conductance. In the two-wire model we have 

Var = ( - j Var Gtot = (^) (Var T, + Var T,) . 

Treating both edges as one-dimensional wires, we map 
the transport problem onto that of a disordered ID 
chain. Transmission eigenvalue statistics in ID disor- 
dered chains is known to be described by the Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation [24 . Using the 
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FIG. 3: Step disorder: edge disorder created by a random 
walk, where the width of the nanoribbon is changed by one 
hexagon at every step. Steps are made with probabihty a/d 
and the maximum deviation of the width is < s hexagons. 
Single vacancies: edge atoms are removed randomly with the 
probability a/d. Extended vacancies: similar to single vacan- 
cies, but also neighboring edge atoms are removed. 




FIG. 4: (color online) Average total conductance, (Ctot) 
(blue solid line), rms of the total conductance, rms Gtot (black 
dashed line), and rms of the spin conductance, rmsCs (red 
solid line), as a function of Ep {Ep = is chosen to corre- 
spond to zero gate voltage). The data were averaged over 1000 
configurations of single vacancies with d = 40a and L = 800a. 
For comparison, the inset shows the same quantities for for 
the singular case of = 0. In this situation, the spin conduc- 
tance and its fluctuations vanish completely. 

full distribution function of resistance [25 , we find that 
the universal maximum value of the root mean square 
(rms) spin conductance rmsGg = VVar Gs ^ 0.4e/47r. 
In order to demonstrate this universality, we investigate 
GNRs of different length L and width W and various 
models of edge disorder (see Fig. [3| . 

First, we focus on dilute disorder, where the average 
distance between scatterers d ^ a. The typical behav- 
ior of charge and spin conductances (average, fluctua- 
tions) is shown in Fig. [4j We first note that over the 
whole energy region, where the edge states are present, 
^ rms Gtot ^ rmsGs, confirming the validity of the two- 
wire model. As the Fermi level is raised by gating or 
doping, the relevant states are extended and feel both 
edges. Then, the assumption of uncorrelated channels 
breaks down, and ^ rms Gtot > rmsGg. 

For an n-type GNR, when the Fermi level is near 
the band edge, the states at Ep are localized and both 
the average conductance and the fluctuations are sup- 



pressed exponentially. Raising E^p, we observe in Fig. [4] 
a crossover to the ballistic regime, where the conductance 
rises up to the quantum limit of conductance 2e^ /h. Cor- 
respondingly, we see a maximum in the conductance fluc- 
tuations before they vanish again in the ballistic regime. 

The average/fluctuations of the conductances of a p- 
doped GNR are different from an n-doped one, but a 
description based on the DMPK equation holds well 
for either case. The scattering strength of impurities 
depends on the overlap of the impurity potential with 
the unperturbed channel wavefunction and therefore on 
Aedge = Aedge(^F)- In the u-dopcd GNR, there is one 
channel whose momentum is a monotonic function of Ep. 
On the other hand, in the p-doped GNR, due to the band 
dispersion (Fig.jl]), there are two channels: One localized 
near the edge, the other extended further into the ribbon, 
but still with a considerable density at the edge. Low- 
ering Ep thus localizes one state even more towards the 
edge, whereas the other state spreads out, making the 
density more uniform. This leads to different functional 
dependences of the localization length on the Fermi en- 
ergy for n- and p-doped ribbons. 

In order to compare n- and p-doped ribbons as well as 
different disorder models, we extract the energy depen- 
dence of the longitudinal (transport) localization length 
^(^f) from exp((/n(G|/|(£^F,I^))) = exp(-2L/0 [M 
(27j, as shown in the inset of Fig. [sja). In Fig. |5|a) we 
show rmsGs as a function of ^/L for all three disorder 
models (see Fig. [3| with different values of d and a wide 
range of ribbon lengths L. The data collapse onto a single 
curve, demonstrating the universality of the spin conduc- 
tance fluctuations (SCF), independent of the particular 
type of edge disorder. Slight deviations from this univer- 
sality can be observed in Fig.jsja), in the ballistic regime 
for the special case of single vacancies. In this case, the 
system reaches the ballistic limit only for high Fermi en- 
ergy values, where the two-wire model breaks down. The 
rms spin conductance of the n-doped GNR agrees very 
well with the results obtained from the DMPK equation. 
For the p-doped ribbon, where there are two conducting 
channels, we see a small increase in the rms conductance, 
presumably due to the crossover to a multi-channel quasi- 
ID wire, where rms G K 0.52 [22j. In Fig.lgb) we con- 
centrate on n-doped graphene for step disorder (upper 
panels of Fig. |3| and show again the universality of the 
SCF with respect to a wide range of parameters charac- 
terizing edge roughness, ribbon length and width. There 
is little dependence on the ribbon width W ^ confirming 
that the observed effect is entirely due to the edges. 

Currently there is not much experimental control over 
the edges of nanoribbons. Considering GNRs with dense 
disorder, d = 0(a), the observed maximum of SCF de- 
creases with increasing disorder density, i.e. decreasing 

as shown in the inset of Fig. jsjb). We observe that 
for d > 5a the SCF are independent of the maximum 
height of the steps. Moreover, we find that the maxi- 
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FIG. 5: (color online) Spin conductance fluctuations: (a) 
rmsCs as a function of ^/L for n- and p-doped graphene: 
step disorder for n-type, d = 20a, s = 3 (black), single vacan- 
cies for n- and p-type, d = 40a (red and blue, respectively) 
and extended vacancies for n-type, d = 30a (green). Inset: 
^/a as a function of Ef for different disorder models (colors 
as in the main panel), (b) rmsGs as a function of ^/L for 
step disorder in n-doped graphene: d = 20a and s = 3 (red; 
orange for W = 92a/ ), = 35a and s = 2 (black), d = 35a 
and s = 6 (blue; violet for W = 92a/ VS), d = 20a and s = 6 
(green). Inset: maximum value of rmsGs as a function oi d/a 
for the step disorder models. In both (a) and (b), the solid 
line corresponds to the DMPK prediction. The data is shown 
for GNR lengths L - 800a (Q), 1000a (□), 1200a (A), 1400a 
(+), and 1600a (x), width W — ?>2/\/?>a unless specified oth- 
erwise. The rms Gs is estimated from 1000 {W = ?>2a/V^) 
and 750 {W — 92a/ v^) disorder configurations. 



mum value of the SCF is retained for (i > 5a. As an 
example, the system depicted in the upper right corner 
of Fig. |3] shows spin conductance ^ 0.4e/47r. The finite 
spin conductance of GNRs predicted above, and thus the 
existence of the edge state magnetism, can be detected 
by measuring charge conductance, e.g. by attaching fer- 
romagnetic leads in a two- or four-probe measurement 
similar to Ref. [5 , with one lead being a zigzag GNR. 

In conclusion, we have discussed the spin transport 
properties of graphene nanoribbons. We have shown that 
an ideal GNR has zero spin conductance but nonzero 
spin Hall conductance. Moreover, only GNRs with im- 
perfect edges exhibit a nonzero spin conductance The 
fluctuations of the spin conductance are universal with a 
maximum rms conductance ^ OAe/Air. Thus, graphene 
nanoribbons can be used as an efficient alternative to 
ferromagnetic leads, paving the way to all-graphene spin- 
tronics devices. 
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